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Abstract 



An analytical and numerical study is presented of transmission of radiation 
through a multi-mode waveguide containing a random medium with a complex 
dielectric constant e = e' + ie" . Depending on the sign of e", the medium 
is absorbing or amplifying. The transmitted intensity decays exponentially 
oc exp(— L/(^) as the waveguide length L oo, regardless of the sign of e" . 
The localization length ^ is computed as a function of the mean free path 
/, the absorption or amplification length |(t|~^, and the number of modes in 
the waveguide N. The method used is an extension of the Fokker-Planck 
approach of Dorokhov, Mello, Pereyra, and Kumar to non-unitary scattering 
matrices. Asymptotically exact results are obtained for » 1 and \a\ » 
1/N'^I. An approximate interpolation formula for all a agrees reasonably well 
with numerical simulations. 

PACS numbers: 78.45. +h, 42.25.Bs, 72.15.Rn, 78.20. Ci 
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I. INTRODUCTION 



Localization of waves in one-dimensional random media has been studied extensively, 
both for optical and for electronic systems.0l An analytical solution for the case of weak 
disorder (mean free path / much greater than the wavelength A) was obtained as early as 
1959 by Gertsenshtein and Vasil'ev.i The transmittance T (being the ratio of transmitted 
and incident intensity) has a log-normal distribution for large lengths L of the system, with 
a mean (InT) = —L/^ characterized by a localization length ^ equal to the mean free path. 

This early work was concerned with the propagation of classical waves, and hence in- 
cluded also the effect of absorption. In the presence of absorption the transmittance decays 
faster, according toSi (InT) = (a — /~^)L, where \a\ is the inverse absorption length [a < 0). 
Absorption is the result of a positive imaginary part e" of the (relative) dielectric constant 
e = e' + is" . For a homogeneous e" one has 



a 



-2k Im Vl + i£" ^ -ke" if \e"\ < 1, (1.1) 



where k is the wavenumber. A negative e" corresponds to amplification by stimulated 
emission of radiation, with inverse amplification length a > 0. Propagation of waves through 
amplifying one- dimensional random media has been studied in Refs. |5|-[T(]|. In the limit 
L ^ oo amplification also leads to a faster decay of the transmittance, according to (InT) = 

A natural extension of these studies is to waveguides which contain more than a single 
propagating mode. Localization in such "quasi-one-dimensional" systems has been stud- 
ied on the basis of a scaling theory,0 a supersymmetric field theory,El or a Fokker-Planck 
equation.lli'0 It is found that the localization length for modes is enhanced by a factor of 
order A^ relative to the single-mode case. These investigations were concerned with quantum 
mechanical, rather than classical waves, and therefore did not include absorption. It is the 
purpose of the presentpaper to extend the Fokker-Planck approach of Dorokhov, Mello, 
Pereyra, and KumarEl'Ei (DMPK) to include the effects on the transmittance of a non-zero 
imaginary part of the dielectric constant. 

According to the general duality relation of Ref. |^, the localization length is an even 
function of a for any A^, 

e(a) = e(-a). (1.2) 

It follows that both absorption and amplification lead to a faster decay of the transmittance 
for large L. For A^ ^ 1 we find that, in good approximation. 

This result becomes exact in the two limits |cr| ^ 1/N'^l and \a\ ^ 1/N'^l. We compare 
with numerical solutions of the Helmholtz equation, and find reasonably good agreement 
over the whole range of a. 

The outline of the paper is as follows. In Sec. II we formulate the scattering problem and 
summarize the duality relation of Ref. ^. In Sec. Ill we derive a Fokker-Planck equation for 
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the transmission and reflection eigenvalues 7^, Tin, n = 1,2,. . .N. These are eigenvalues 
of the matrix products tt^ and rr^ , respectively, where t and r are the transmission and 
reflection matrices of the waveguide. For a = the Fokker-Planck equation is the DMPK 
equation.lli'0 A reduced Fokker-Planck equation, containing only the TZnS, was previously 



obtained and studied in Ref. 15. To obtain the localization length one needs to include 



also the T^'s, which are no longer related to the 7^„'s when a ^ 0. We find that a closed 
Fokker-Planck equation containing TZ^s and 7^'s exists only for N = 1. If > 1 there 
appears an additional set of "slow variables," consisting of eigenvectors of rr^ in a basis 
where tt^ is diagonal. (These new variables do not appear when a = 0, because then rr^ 
and tt^ commute.) Because of these additional relevant variables we have not been able 
to make as much progress in the solution of the Fokker-Planck equation for a 7^ as one 
can for a = 0.0 In Sect. IV we show that a closed evolution equation for (InT) can be 
obtained if |cr| ^ 1/N'^l, which leads to the second term in Eq. (|1.3| ). (This term could 
also have been obtained from the incoherent radiative transfer theory for o" < 0, but not 
for cr > 0.) To contrast the multi-mode and single- mode cases, we also briefly discuss in 
Sec. IV the derivation of the localization length for = 1. (Our A^ = 1 results were given 
without derivation in Ref. ^.) Finally, in Sec. V we compare the analytical results for the 
multi-mode case with numerical simulations. 



II. FORMULATION OF THE SCATTERING PROBLEM 

We consider a random medium of length L with a spatially fluctuating dielectric constant 
e = e' + is" , embedded in an A^-mode waveguide with e = 1. The scattering matrix S is 
a 2N X 2N matrix relating incoming and outgoing modes at some frequency uj. It has the 
block structure 

S=(;' t'). (2-1) 

where t, t' are the transmission matrices and r, r' the reflection matrices. We introduce the 
sets of transmission and reflection eigenvalues {T^}, {T^}, {"^n}, {"^n}, being the eigenvalues 
of, respectively, tt^, t't'^, rr\ _ Total transmittances and reflectances are defined as 

T = N-^i:iU\ R = N-^TTrr\ (2.2a) 
r = Trt't't, R' = N-^ Tir'r'l (2.2b) 

Here T and R' are the transmitted and reflected intensity divided by the incident intensity 
from the left. Similarly, T' and R correspond to incident intensity from the right. By taking 
the trace in Eq. ( |2.2| ) we are assuming diffuse illumination, i.e. that the incident intensity 
is equally distributed over the A^ modes. Two systems which differ only in the sign of e"{r) 
are called dual. Scattering matrices of dual systems are related byi 

S{e")S^{-e") = 1. (2.3) 

This duality relation takes the place of the unitarity constraint when e" 7^ 0. 
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An optical system usually possesses time-reversal symmetry, as a result of which 
Sie")S*{-e") = 1. Combining this relation with Eq. (U), we find that ^ = is a 
symmetric matrix. Hence Tn = and T = T' . (The refiectances R and R' may differ.) The 
case of broken time- reversal symmetry might also be physically relevant ,0 and is included 
here for completeness. In the absence of time-reversal symmetry S is an arbitrary complex 
matrix. 

The duality relation ( |2.3| ) has consequences for the reflection and transmission eigenvalues 
of two dual systems.! If = 1 the relation 

T{e")/R{e") = r{~e")/R'{-e") (2.4) 

holds for all L. If > 1 we have two relations for L —>■ oo, 

lim 7^„(£") = lim n-\-e"), (2.5) 

L— >oo L— >oo 

lim L-Hn%{6") = lim L-Hn%{-e"). (2.6) 

L—*oo L— >oo 

The transmittance T = N~^J2n%i is dominated by the largest transmission eigenvalue, 
hence 

lim L~MnT(£") = lim L-HnT{-6"). (2.7) 

>oo L—*oo 

In other words, two dual systems have the same localization length, as stated in Eq. ( |1.2| ). 



III. FOKKER-PLANCK EQUATION 

We derive a Fokker-Planck equation for the evolution of the distribution of scattering 
matrices with increasing length L of the waveguide. In the absence of gain or loss (cr = 0)^he 
evolution equation is known as the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation.ll^ll3 
Original derivations of this equation relied on the unitarity of the scattering matrix, making 
use of the invariant measure on the unitary group and the polar decomposition of a unitary 
matrix. These derivations cannot readily be generalized to the case a 7^ 0, in particular 
because the scattering matrix no longer admits a polar decomposition. (This means that 
the matrix products rr'^ and tt^ do not commute.) The alternative derivation of the DMPK 
equation of Ref. |1^ does not use the polar decomposition and is suitable for our purpose. 

Without loss of generality we can write the transmission and reflection submatrices of 
the scattering matrix as follows, 

^_/r' t'\_fUVR^W U'VT^Z \ 

^~[t r)~[vVTW' -V'VKZ'J- ^"^-^^ 

Here U, U', V, V, W, W, Z, Z' areNxN unitary matrices, while R, R', T, T' are diagonal 
matrices whose elements are the reflection and transmission eigenvalues {7^„}, {7?.^}, {%i}, 
{X}. For a = 0, the unitarity constraint SS^ = 1 implies U = U' , V = V , W = W, Z = Z', 
and R = R' = 1 — T = 1 — T'. Eq. ( PTT| ) then constitutes the polar decomposition of the 
scattering matrix. In this case one can derive a Fokker-Planck equation for the evolution 
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of only transmission or only reflection eigenvalues. If o" 7^ 0, the Fokker-Planck equation 
contains both the transmission and reflection eigenvalues, as well as elements of the matrix 
Q = V'^V relating eigenvectors of tt^ and rr^. The only constraint on the scattering matrix 
if cr 7^ is imposed by time- reversal symmetry, which requires S = S"^, hence W = U"^, 
Z = V^, W = U'^, Z' = V'^, T = T'. 

The Fokker-Planck equation describes the evolution of slow variables after the elimination 
of fast variables. In our problem fast variables vary on the scale of the wavelength A, while 
slow variables vary on the scale of the mean free path / or the amplification length |cr|~^. 
We assume that both / and are much greater than A. (This requires \e"\ <C 1.) The 
slow variables include {TZn}, {^} and elements of Q = V^V. We denote this set of slow 
variables collectively by {$n}- Each $j is incremented by if a thin slice of length SL 
{X <^ SL I) is added to the waveguide of length L. The increments are of order {SL/lY^"^ 
and can be calculated perturbatively. We specify an appropriate statistical ensemble for the 
scattering matrix SS of the thin slice and compute moments of 5$^. The first two moments 
are of order SL/l, 

{6<l>,) =ai6L/l + 0{6L/lf/^, (3.2a) 
(5$i5$j)= aij 5L/1 + 0{5L/lf/^. (3.2b) 

Higher moments have no term of order 5L/1. According to the general theory of Brownian 
motion,El the Fokker-Planck equation for the joint probability distribution P({$„}, L) reads 

dP d 1 (9^ 

= -Y w-<^iP + -Y w^^ciiiP- (3.3) 



The average (■ ■ ■) in Eq. ( p^ ) is defined by the statistics of 5S. We specify this statistics 
using simplifying features of the waveguide geometry (length ^ width), which justify the 
equivalent-channel or isotropy approximationJli'S We assume that amplification or absorp- 
tion in the thin slice is independent of the scattering channel. This entails the relation 

5S5S^ = l + a5L, (3.4) 

where a is a modal and spatial average of the inverse amplification length a. If e" is spatially 
constant, one has 

-=-|El-(l--n/^' + ^^")^^ (3.5) 
n=l 

where ujn is the cutoff frequency of mode n. For iV — >■ 00, the sum over modes can be 
replaced by an integral. The result depends on the dimensionality of the waveguide, 

a = —2ke", for a 3D waveguide, (3.6a) 
a = —{TT/2)ke", for a 2D waveguide, (3.6b) 

where we have used that \e"\ <^ 1. 

Eq. (|3.4| ) ensures the existence of a polar decomposition for 6S, 



5 



6S 



U0VM.W0 
VoVSTWo 



UoVSTZo 
-VoVSRZo 



(3.7) 



with 6T + 5R = 1 + a6L. Note that a polar decomposition for 6S does not imply a 
polar decomposition for S, because the special block structure of Eq. ( |3.7| ) is lost upon 
composition of scattering matrices. We make the isotropy assumption that the matrices Uq, 
Vq, Wq, Zq are uniformly distributed in the unitary group. In the presence of time-reversal 
symmetry one has Wq = Uq and Z^ = V^. In the absence of time- reversal symmetry all four 
unitary matrices are independent. The diagonal matrices (5R and 6T may have arbitrary 
distributions. We specify the first moments, 



(Tr5R) = NSL/l, 

{TtST) = N + N{-f -l)6L/l, 

where we have defined 'y = al. The mean free path / in Eq. 
path /ti- of radiative transfer theory by@ 

/ = (4/3)/tr, for a 3D waveguide, 
I = (7r/2)/ti., for a 2D waveguide. 



(3.8a) 
(3.8b) 

is related to the mean free 

(3.9a) 
(3.9b) 



This completes the specification of the statistical ensemble for 6S. 

We need the increments A7^„, A7^ of reflection and transmission eigenvalues to first 
order in 6L/1, 

Ai?« + AR^l + E (3.10a) 



TZ 



AT(i)AT(i) 

nm mn 

T - T 



(3.10b) 



The matrices of perturbation AR^^\ AR^'^\ AT^^\ AT^^-* are expressed through unitary 
matrices Q = V^V, U = Z'Uq, W = WqV and diagonal matrices T, R, 6T, (5R, 

Ai?(^) = \^/RUV6RW{l - R) + H.c] , (3.11a) 

(3.11b) 
(3.11c) 



- '^W\l - 6T)WR + VWV6RWVkuV6RW{1 - R) + H.c. 
AT(i) = - IqVkuVSRWQ^T + H.c 



AT(2) = Q\^uVmWQ^TQW^V6RWVKQ^ 

- [^QW\l - 6T)WQ^T - QVrU\^wVkuV6RWQ^T + H.c] . (3.11d) 

(The abbreviation H.c. stands for Hermitian conjugate.) The moments (|3.2|) are computed 
by first averaging over the unitary matrices Uq, Wq and then averaging over 5R and 6T 
using Eq. ( p.8|) . Averages over unitary matrices follow from 

1 



N 



' ^nm ^kl ) 



{^nk^mk^pl^ql) ~ _/\r^_/\r _|_ 1)*^^' {^rip^mq + 5nq5mp) 



(3.12a) 
(3.12b) 
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Without time- reversal symmetry averages over Uo and Wo are independent. With time- 
reversal symmetry we have Wq = Uq so that only a single average remains. The results 
are 



(l/5L){5nn) = 1 + 2(7- l)7^, 
1 



+ 



With time-reversal symmetry 
N + 

7^„(l-7^J2 + 7^^(l_7^J2 



N + 



{i/SL){snjnm) 



45. 



N + l 



nJi-TZn] 



{l/5L){5%) = T„(7 - 1) + j^{Ann + E -^"^^"^ ^ 



7-7 



7+7 



{l/6L){S7J7,n) 
{l/5L){57jnm) 



N + 



-7 71 (1 - TZ )\0 1^ 



N+l 

Without time-reversal symmetry 



{1/5L) (57^„) = 1 + 2(7 - 1)7^„ + -f E ^ 



1 7^n(l-7^^)^' + 7^^(l-7^„.)^ 



(//5L)(57^„57^„) 



2(5„ 



-7^„(l - Tin] 



(//5L)(5T„) = T„(7 - 1) + |(a.„ + E -^"^^"^ ^ ^"^^ 



7-7 



{l/6L){d7J7J 
{l/6L){67jn^) 



25, 



N 
2 



-'n ^nni 



(3.13a) 
(3.13b) 

(3.13c) 

(3.13d) 
(3.13e) 



(3.14a) 

(3.14b) 
(3.14c) 

(3.14d) 
(3.14e) 



We have abbreviated Anm = (QR-Q"^)nm and Fnm = \{Qy/R.Q'^)mn\'^- 

The moments of 6TZn contain only the set of reflection eigenvalues {TZn}, so that from 
Eq. (3^) we can immediately write down a Fokker-Planck equation for the distribution of 
the T^n's. In terms of variables /i„ = l/(7^„— 1) G {—oo, —1) U (0, oo) it reads@ 



d 



N 



d 



E ^ — /""(I + 



dP 



1 



(3.15) 
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where the symmetry index (3 = 1(2) corresponds to the case of unbroken (broken) time- 
reversal symmetry. The evolution of the reflection eigenvalues is independent of the trans- 
mission eigenvalues — but not vice versa. The evolution of the 7^'s depends on the TZnS, 
and in addition on the slow variables contained in the unitary matrix Q. To obtain a 
closed Fokker-Planck equation we also need to compute increments and moments of Q. The 
resulting expressions are lengthy and will not be written down here. 

In the single-mode case (A^ = 1) this complication does not arise, because Q = e'^'^ 
drops out of the scalars A and F. The single transmission and reflection eigenvalues T, TZ 
coincide with the transmittance and reflectance T, R defined by Eq. ( [2.21 ). The resulting 
Fokker-Planck equation i^ 



—n-, -1 + R)P+ ^T'-RP - 2g^TR(l - R)R (3.16) 



In the case of absorption (7 < 0), Eq. ( p.l6|) is equivalent to the moments equations of 
Ref. |. 

IV. LOCALIZATION LENGTH 

The limit L — * 00 of the distribution of the reflection eigenvalues follows directly from 
Eq. (|3.15| ), by equating the left-hand-side to zero. The resulting distribution Pqo is that of 
the Laguerre ensemble of random matrix theory,lii 

P{{f^n}) oc n - ^^^fl[^M-liPN+2-(3)^,k]. (4.1) 

i<j k 

The distribution looks the same for both signs of 7, but the support (and the normalization 
constant) is different: /i„ > for 7 > 0, and /i„ < —1 for 7 < 0. To determine the 
localization length we need the distribution of the transmission eigenvalues in the large-L 
limit. We consider the cases = 1 and ^ 1. 



A. Single-mode waveguide 

We compute the distribution P{T, L) of the transmission probability through a single- 
mode waveguide in the limit L 00. In the case of absorption (7 < 0) this calculation was 
done by Rammal and Doucot,i and by Freilikher, Pustilnik, and Yurkevich.@ We generalize 
their results to the case of amplification (7 > 0). The two cases are essentially different 
because, while the mean value of R is finite in the case of absorption, 

(i?)oo = 1 - 276-'^ Ei(27), for 7 < 0, (4.2) 

it diverges in the case of amplification. The mean value of In R is finite in both cases, 

/mm _/C + ln27-e2^Ei(-27), for 7 > 0, 

^^^~~\-C-ln(-27) + e-2^Ei(27), for 7 < 0. ^ 
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Here C is Euler's constant and Ei(x) = J^^dte^/t is the exponential integral. The relation 

(lni?(7))oo = -(lni?(-7))oo (4.4) 



holds, in accordance with the duality relation ( |2.5| ). 

We now show that the asymptotic L oo distribution of T is log-normal, with mean 
and variance of In T given by 

(In T) = -(1 + |7|) L/l + 2c(7) + 0(//L), (4.5a) 

/' N _ / 0' for 7 < 0, , , . 

^17) - |c^lj^27-e2^Ei(-27), for 7 > 0, ^ 



var In T 



2 + 4|7|e2|^lEi(-2|7|) L/1 + 0{1). 



(4.6) 



The constant 0(7) ^ -27 In 7 if < 7 < 1. Note that var In T < (Inr)^ for L/l > 1. The 
localization length ^ = Z (1 + |7|)~^ is independent of the sign of 7, in accordance with the 
duality relation (|1.2|) . 

These results are easy to establish for the case of absorption, when Eq. (|3.16| ) implies 
the evolution equation^ll 



l-^ var In T 
oL 



2{R), 



for 7 < 0. 



Making use of the initial condition T 
one readily obtains Eqs. 



1 for L ^ and the asymptotic value 



(4.7) 
_ oi{R), 

5D and (g^) for 7 < 0. 
In the case of amplification, the evolution equations ( [4. 71) hold only for lengths L smaller 
than Lc — /c(7)/|7|. For L ^ Lc stimulated emission enhances transmission through the 
waveguide. On larger length scales stimulated emission reduces transmission. Technically, 
the evolution equations ( |4.71 ) break down for L — > 00 because the integration by parts of the 
Fokker-Planck equation produces a non-zero boundary term if L > Lc. To extend Eqs. ( |4.5|) 
and ( |4.6| ) to the case 7 > we use the duality relation ( p.4| ). It implies that for = 1 the 
distribution of the ratio T/R is an even function of 7. Eq. ( [4.5|) for 7 > follows directly 
from the equality 

(lnr(7)/i?(7)) = (lnT(-7)/i?(-7)), (4.8) 

which holds for all L, plus Eq. ([4.4| ), which holds for L 00. The constant 0(7) for 7 > 
equals (lni?(7))oo and is substituted from Eq. ( ^4.3| ). The duality of T{j)/R{j) also implies 
Eq. (|4.6|) for the variance, provided the covariance ((InTlni?)) = (lnTln/2) — (lnT)(lni?) 
remains finite as L — > cxd. We have checked this directly from the Fokker-Planck equation 
( P.16| ), and found the finite large-L limit 



((lnTlni?))oo = -2e2^Ei(-27)c(7) -c(7)' 



27 / dfx ln^(l + /i) - InV , for 7 > 



(4.9) 
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B. Multi-mode waveguide 



We next consider a waveguide with ^ 1 modes. We compute the localization length 
^ = — limi^oo -^"^(InT) in the case of absorption, and include the case of amplification 
invoking duality. For absorption the average reflectance {R) = A^"^(Eyfc(l + l/z^fc)) 

remains 

finite as L ^ oo. The large-L limit (-R)oo follows from the distribution ( |4.1| ), using known 
formulas for the eigenvalue density in the Laguerre ensemble.@ For |7|A^^ ^ 1 the result is 



(R)^ = 1 + |7l - \/|7l(2 + |7l) + 0{1/N), 7<0. (4.10) 

The evolution of transmission eigenvalues is governed by the Fokker-Planck equation 
with coefficients given by ( ^.2| ), ( ^.131 ), and ( |3.14[ ). Each 7^ has its own localization 



length C,n = — liniL^oo In 7^. We order the ^^'s from large to small, > ^2 > • • • > ^iv- 
This implies that for L — 00 the separation of the 7^'s becomes exponentially large, Ti 3> 
T2 ^ . . . ^ T]\f. Hence we may approximate 

%i + %n f-l, for n > m, \ 

(4.11a) 



Tn-Tm 1 1, forn < m, 

"Tn-^mm ~\~ Tm-^nn ( Ann, fo^ ^ ^ lib) 

The Fokker-Plank equation ( |3.3| ) simplifies considerably and leads to the following equation 
for the largest transmission eigenvalue: 

^^^^"^^^ " I -1 - l7l + {R) - M^n). for (3 = 2. ^^''^^ 

For |7|A^^ ^ 1 we may substitute Eq. (|4.10| ) for (R) and omit the terms with (An) and 
(Fii). The resulting localization length is given by 



= VM2TM) + 0{l/N). (4.13) 

Because of duality, Eq. ([4.13|) holds regardless of the sign of 7. It agrees with radiative 
transfer theory for 7 < 0, but not for 7 > 0. Indeed, the exponential decay of the transmitted 
intensity in the case of amplification is an interference effect, which is not contained in the 
theory of radiative transfer. 

Eq. ( |4.13| ) is asymptotically exact for I7I 3> 1/N'^. For smaller I7I we cannot compute ^ 
rigorously because the distribution of the matrices A and F is not known. An interpolative 
formula for all 7 can be obtained by substituting for (An) and (Fu) in Eq. ( [4.12] ) their 
L — s> 00 limits when 7 = 0, which are (An) = (Fn) = 1. In this way, we arrive at the 
localization length 



pN+2-p 



+ VM2 + H) 



(4.14) 



which interpolates between the knownSliil value of ^ for 7 = and Eq. ( [4.13|) for I7I 3> 

1/N\ 
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The localization length ^ is the largest of the eigenvalue-dependent localization lengths 
^n- What about the other ^„'s? For 7 = it is knownEl'Sill that the inverse localization 
lengths are equally spaced, and satisfy the sum rule J2n VCn = N/l. We have not succeeded 
in deriving the spacings for 7 7^ 0, but we have been able to derive the sum rule from the 
Fokker-Planck equation (by computing the L-dependence of {J2n^^%i))- The result is exact 
and reads 

N 

ij:c' = a+b\)N. (4.15) 

n=l 



V. NUMERICAL RESULTS 

To test the analytical predictions on a model system, we have numerically solved a 
discretized version of the Helmholtz equation, 

+ k'^e{r)]E{r) = 0, (5.1) 

on a two-dimensional square lattice (lattice constant d, length L, width W). The real part e' 
of the dielectric constant was chosen randomly from site to site with a uniform distribution 
between 1 ± 6e. The imaginary part e" was the same at all sites. The scattering matrix was 
computed using the recursive Green function technique.il 

The parameter a is obtained from the analytical solution of the discretized Helmholtz 
equation in the absence of disorder {6e = 0). The complex longitudinal wavenumber kn of 
transverse mode n then satisfies the dispersion relation 

cos{knd) + cos{mTd/W) = 2- l{kdf{l + is"), (5.2) 

which determines a according to a = —2N~^ Im J2n ^n- Simulations with e" = were used 
to obtain I, either from the large-L relationlii 



- lim L-i(lnT) = [1{N + (5.3) 



or from the large- relatione 



lim(T) = (l + L//)-\ (5.4) 

N^QO 



The parameters chosen were W = 25 d, k = 1.22 d ^, corresponding to = 10, / = 29.6 d 
from Eq. (|5.3|) and / = 26.1 d from Eq. ( |5.4|) . The localization length was computed as a 



function of a from the L-dependence of InT up to 40 /, averaged over 150 realizations of the 
disorder. Results are shown in Fig. |^. The localization length is the same for absorption 
and amplification, within the numerical accuracy. Comparison with the analytical result 
( [4.14| ) for /5 = 1 is plotted for the two values of the mean free path. The agreement 



is quite reasonable, given the approximate nature of Eq. ( [4.14|) in the regime |7|A^^ — 1 
(corresponding to |(T|(i ~ 10~^) 
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Localization length ^ = — lim^^^oo ^^^(InT) of a disordered waveguide (A^ = 10) 
versus the modal average a of the inverse absorption or amplification length. Data points are a 
numerical solution of the discretized (lattice constant d) two-dimensional Helmholtz equation for 
the case of absorption (squares) or amplification (circles). The curves are the analytical prediction 
(P^ ) in the case /3 = 1 (unbroken time-reversal symmetry) for / = 29.6 d [solid curve, determined 
from Eq. (5^)] and for / = 26.1 (i [dashed curve, determined from Eq. (|5.4| )]. The inset shows the 
same data on a linear, rather than logarithmic, scale. 



In conclusion, we have shown how absorption or amplification can be incorporated into 
the Dorokhov-Mello-Pereyra-Kumar equation for the transmission through a multi-mode 
waveguide. The technical difficulty of the multi-mode case is that the Fokker-Planck equa- 
tion for the transmission eigenvalues 7^ depends not just on the transmission and reflection 
eigenvalues 7^, 7^^, but also on the eigenvectors of the matrices tt'^ and rr^'. We could com- 
pute the localization length in the two regimes, I7I 3> and I7I -C and have given 
an interpolation formula for the intermediate regime. An exact solution for all 7 remains an 
unsolved problem. 
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